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We develop an approach describing nonlinear-optical processes in the strong-held domain character¬ 
ized by the nonperturbative held-with-matter interaction. The polarization of an isolated atom in 
the external held calculated via the numerical solution of the time-dependent Schrodinger equation 
agrees with our analytical hndings. For the practically important case of one strong laser held and 
several weaker helds we derive and analytically solve propagation equations describing high-order 
(HO) wave-mixing, HO parametric amplihcation and HO stimulated scattering. These processes 
provide a way of efficient coherent XUV generation. Some properties of HO processes are new in 
nonlinear optics: essentially complex values of the coefficients in the propagation equations, the 
super-exponential (hyperbolic) growing solutions, etc. Finally, we suggest conditions for the prac¬ 
tical realization of these processes and discuss published numerical and experimental results where 
such processes could have been observed. 
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The nonlinear optics usually deals with the laser field 
which perturbatively interacts with matter. In this case 
the matter response is described with nonlinear suscepti¬ 
bilities • This approach allows investigating numer¬ 
ous nonlinear optical effects involving few photons [1]. 
However, in the case of intense laser field the perturba¬ 
tion approach fails. A number of multiphoton process in¬ 
volving electronic dynamics (ionization, electronic rescat¬ 
tering, etc) are successfully described within nonpertur¬ 
bative approaches (like Keldysh approximation [5]), but 
the only fully optical process which is well-understood in 
this case is high-order harmonic generation (HHG). The 
study of the other nonlinear optical processes in the non¬ 
perturbative regime is limited, among other factors, by 
poorly developed theoretical methods of their descrip¬ 
tion. For instance experimental [3] observation of the 
exponential growth of the high harmonic signal has led 
to an active discussion gHS]. The feasibility of high-order 
optical processes was already discussed in the early stud¬ 
ies [7HTQ] but the theoretical methods were based on us¬ 
ing nonlinear susceptibilities and thus were limited to 
processes of relatively low order. More recently the XUV 
amplification was obtained in simulations m based on 
direct numerical integration of the propagation equation 
with the nonlinear polarization calculated via numeri¬ 
cal solution of the time-dependent Schrodinger equation 
(TDSE). 

In this paper we suggest an approach to describe non¬ 
linear optical processes in the presence of a given strong 
laser field denoted as Eq. Processes involving other helds 
(denoted as Ei) are described with the susceptibilities, 
nonperturbatively induced by the pump held. Irrespec¬ 
tive of the nature of nonlinearity, the microscopic polar¬ 
ization P{t) is a functional of the external held: P(t) = 
<&[!?] = + El], The functional can be expanded 

in the Taylor series |12j (applications in fundamentals 
of quantum mechanics |I3) and density functional the¬ 
ory [14j are well-known) using so called functional deriva¬ 
tives of the functional $ over the function E. This ex¬ 


pansion includes integration and thus differs from the 
traditional Taylor expansion of a function. The details 
of the calculation are presented in Appendix I. Consider¬ 
ing the helds Agp = £o,i exp(—jwoxt+ kopr-l-(/:o,i) + c.c. 
with the real amplitudes Sop, we expand the functional 
derivatives, as well as the functional in the Fourier 

series. Finally we hnd that 

exp{iq{—uJot -b kor + ipo) + im{—uJit + kir + ipi)} + c.c. 

( 1 ) 

where the induced susceptibilities are 

+ 00 

ltM(£-o,fi,a;o,cci)= ^ 4"^’^-)(£o,cco,a;i)|£iP (2) 

i=o.2,-- 

Thus, the microscopic response is a sum of waves to 
which q photons from the one held and |m| -I- j photons 
from the other one contribute. For the case Wq = Wi, 
ko yf ki similar equation was found in m- Note that 
up to here we never supposed that the held £i is small. 
Terms with j = 2,4,... in equation ([^ can be under¬ 
stood as those describing processes where one, two, etc. 
photons from the held Ei were absorbed and the same 
number of photons of this held were emitted. 

When £i « Sq the hrst term in equation ([^ domi¬ 
nates, so we have k^\£o,£i,u}o,^i) ~ k^’^\£o,ujo,ujA 
As does not depend on £i, in the expansion (IH 

the dependence on the held £i remains only in the term 
(fi)!"*!. This makes this expansion for the polarization 
similar to the one appearing in the usual (’’perturba¬ 
tive”) nonlinear optics [1]. Below for brevity we write 

Am) _ (m.O) 

/\q — tSjq 

The remarkable practical importance of equation 0 
is that even in the strong held domain one can write and 
sometimes solve analytically (see below) the propagation 
equation with the right side given with equation 0 . So 
there is no need to calculate the non-perturbative re- 
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Figure 1: (color online) Numerical spectra of the the polar¬ 
ization induced by the weak field E\ in presence of the intense 
field Eq for different strengths of the weak field, see text for 
more details. Arrows show the spectral components generated 
due to different induced susceptibilities Kq"*'. The spectra are 
normalized with the intensity of the 13-th harmonic generated 
in absence of the weak field. 


sponse simultaneously with the propagation equation so¬ 
lution. Instead one should find the polarization using 
an appropriate theoretical approach (this can be numeri¬ 
cal TDSE solution, strong-held approximation for an iso¬ 
lated atom response, PIC simulations for plasma, etc), 
calculate its spectrum P(a;) and then hnd k from this 
spectrum using equation Q; alternatively, one can di¬ 
rectly calculate n from equations (12)-(14) (see Appendix 

I)- 

To check our analytical hndings we compare the found 
properties of the polarization with numerical results. We 
solve numerically the 3D TDSE in the single-active elec¬ 
tron approximation for a model Xe atom in an exter¬ 
nal laser held, the details of the TDSE solution are pre¬ 
sented in [TB]. For the given strong held Eq we cal¬ 
culate numerically the polarization P in the presence 
of the weak held Ei and the polarization in the 
absence of this held. Figure presents the spectrum 
P(uj,Ei) = P{uj,Ei) — p(°)(w) for the given £q and dif¬ 
ferent Si. The strong held wavelength is 800 nm and 
its intensity is 5 x 10^^ W/cm^. The weak held has low 
frequency uji = |wo and low intensity, which differs in 
different calculations. 


We can see that in agreement with the selection rule 
for a central-inverse medium, only processes with an odd 
number of photons <7 -I- \m\ + j contribute to the pro¬ 
cess (for \m\ < 2 and j = 0 this selection rule was 
demonstrated already in the early experiment |17| ). The 
strongest contributions are the ’’doublet” around the 
even harmonic of the pump held g = 12 due to the 
|m| = 1 , j = 0 processes , and the ’’triplet” around the 
odd one g = 13 due to the \m\ + j = 2 processes. 

Figure shows the absolute value and the phase of 



Figure 2: (color online) Absolute values and phases of spectral 
components of the response shown in Fig as a function of 
the weak field strength. Several components with different 
frequencies are shown, the induced susceptibility responsible 
for the generation of every component is presented in square 
brackets. Dotted lines show the power approximations. 


the response P{uj,Si) as a function of the held Si. We 
can see that in agreement with our analytical hndings 
the arguments are approximately constant and the abso¬ 
lute values increase approximately as , with low¬ 

est possible j. However, the deviation from the approx¬ 
imations increase with Si because terms with higher j 
should be taken into account. The role of terms with 
J 7 ^ 0 can be seen in this hgure from the component with 
uj = Ihwo (red line with circles) generated due to 
For Si = 0.1£o we can see that |P(13wo)| ~ |P^°^(13ciJo)|, 
so this contribution is not negligible for such Si. Below 
we shall deal with lower helds Si , thus studying only the 
processes with j = 0 . 

In Appendix II we derive some properties of the in¬ 
duced susceptibilities, in particular, their permutation 
symmetry. 

Below we study several high-order nonlinear processes 
illustrated with Figure 

Wave-mixing is the process where g photons from the 
intense field and m photons from the weaker field are con¬ 
verted in one W 2 = + rnuji photon. From equation 

([^ one can see that the intensity of the generated field S 2 
is proportional to the m-th power of the intensity of the 
field Si. This explains the experimental findings of [18] 
where for XUV generation this behavior was observed up 
to TO = 6 . Another example is generation of THz radi¬ 
ation using a frequency-tunable half harmonic of a fem¬ 
tosecond pulse [19]. This process can be understood as 
difference frequency generation (one fundamental photon 
plus TO = — 2 half harmonic photons). The experimen- 
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Figure 3: (color online) The schematic of several high-order 
processes involving many photons from the pump field with 
frequency wo and few photons from the weaker fields with 
frequencies uji and (jJ 2 ■ 


(a) 


(b) 


(c) 




Figure 4: (color online) Intensities of the fields Ei (dotted red 
line) and E 2 (solid blue line) as a function of the propagation 
length for: two-photon parametric amplification (a), three- 
photon parametric amplification for different initial intensities 
of the fields (b and c). 


tally found dependence of the THz yield on the weak half 
harmonic intensity is quadratic, as it should be for the 
process with \m\ = 2 . 

From equation 0 one can see that the detuning 
from the phase-matching for a process of mixing of 
q photons from one beam and m photons from the 
other is, irrespective of the intensity ratio: = 

qkg + mki — k 2 where ^2 is the wave-vector at the fre¬ 
quency W 2 . In m the phase-matching for the non- 
collinear HHG (wi = ojq, ki 7 ^ kg) was considered. 
Here we assume uji 7 ^ wg, co-directed ki and kg, plasma 
and/or capillary contribution to the dispersion domi¬ 
nates, the plasma frequency ujpi << wg, wi. Then 

Afc,,™ = + ■ Omitting the last 

term we find that for w = — ^ the detuning is zero irre¬ 
spective of the plasma frequency, thus irrespective of its 
density. This is very important because the medium ion¬ 
ization always accompanies processes in strong laser field 
and thus the plasma density is time-varying. 

The possibility of the phase-matching optimization in 
the difference-frequency mixing in plasma, involving few 
photons from the two waves was first shown in nnun] 
and further studied in [22]. For the case wg = 2a;i, 
q = 6, m = —3 the phase-matching was experimentally 
demonstrated in |23| . However, the case of wg = 2uJi is 
hardly perspective for solving phase-matching problem 
for really high harmonics because along with the opti¬ 
mal (m = —g/ 2 ) polarization wave, many other polar¬ 
ization waves are generated. This is not the case when 
only few waves with different (small) m are generated. 
For |m| << Igl it should be wi << wg, that is why in 
Fig. [3 and [3 we present the numerical results for this 
case. To achieve phase-matched generation of XUV us¬ 
ing given number q and lowest |m|, for the case of even 
q one should use frequency wi = ujo/q and m = — 1 pro¬ 
cess, and for odd q one should use frequency uji = 2u!o/q 
and m = —2 process, see FigjU 


Parametric amplification and generation. The pro¬ 
cess of the parametric amplification is described with the 
same equations as the difference-frequency generation, 
the difference is in the input conditions: the former as¬ 
sumes that initially there is one intense pump field and 
both generated helds are initially weak, while the latter 
assumes two initial intense fields [T|. 

(i) generation of two photons. Let us consider the pro¬ 
cess in which q photons from the initial field are converted 
into one wi photon and one UJ 2 photon, see Fig. 0 - In 
this case q is an even number (note that for q = 2 this is 
a well-known process of four-wave mixing; however, for 
higher q the analogy with the four-wave mixing is hardly 
helpful). For m = — 1 we have from equation Q: 

= <7Wg - a;2,i)x .gx 

exp{iq{koz - wgt -b (po)}E^ .^ 

Substituting this polarization into the propagation equa¬ 
tion we have: 

881 ^ 2 !dz = *27rA:y2KF)(u;i_2)'?2.i expjfA/cz -b iqipo} (4) 

where in the considered case is Afc = qko — ki — k 2 . 

These equations are similar to those describing the 
parametric amplification in the perturbative nonlinear 
optics (see [Tj, part 9.1). However, an important dif¬ 
ference from this case is that the susceptibility is 
complex, and its phase is not negligible. For instance, 
the phase of is well-understood within the recolli¬ 
sion picture (see [23] and references therein). The details 
of the solution of the propagation equations Q are pre¬ 
sented in Appendix HI. We obtain the solution in which 
the exponential growth dominates after a certain propa¬ 
gation distance. 

Our numerical calculations show that within the 
plateau region is almost independent on q. For 

the conditions of Fig. we find that for q = 12 = 

4.5 X 10~^^ngas CGS units, where Ugas is the gas density 
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(in cm ^). According to equation (35) (see Appendix III) 


this gives the exponential growth of the intensity with the 
increment g = 1.35 cm“^ for the atmospheric pressure. 
Fig. ^a) shows the fields’ intensities calculated for these 
conditions, under zero incident intensity of the first field. 
One can see that initially {z <1 cm) the intensity of the 
first field growths quadratically, so the solution describes 
difference frequency generation, and for z > 2 cm it grows 
exponentially. Note that the solution in the transition re¬ 
gion can be less smooth depending on the initial phases 
of the fields and phases of the induced susceptibilities. 

The found values of increment g show that the high- 
order parametric process can hardly be observed in 
atomic gas jets, but it can be observed in larger targets 
(capillaries and cells) especially under high gas pressure 
(pressures up to tens of atmospheres were used in recent 
experiments |25j 1. Moreover, this process can be very ef¬ 
ficient in parametric generator. Such parametric genera¬ 
tor can be only singly-resonant one because of extremely 
low efficiency of optical elements in XUV. Note that the 
intracavity HHG was demonstrated experimentally [26] 
and now it is actively used. 

(ii) generation of three photons. Another parametric 
process, which was not considered in the perturbative 
nonlinear optics, is the one where q photons from the 
pump field are split into more than two photons. Here we 
shall study the process where these are two low-frequency 
photons (wi) and one XUV photon {uj 2 ). Here q should 
be odd, so under wi << uiq the generated XUV frequency 
0 J 2 is close to the q-th harmonic frequency. Thus, the nat¬ 
urally broadened HH line can provide an effective seeding 
for this process, and this determines its practical impor¬ 
tance. In Appendix IV we obtain the following propaga¬ 
tion equations describing this process: 

dSijdz = AiSlEtf exp{lAA: 2 : -|- iqipo} 

882 !dz = A 2 (f*)^ exp{iAA:z -I- iqtpo} 


where lS.k = gfco — 2/ci — ^ 2 . We solve it analytically 
for Afc = 0. The main property of the found analytical 
solution is the hyperbolic growth after a certain propaga¬ 
tion distance: £’ 1 ^ 2 ( 2 ) c*: 1/(1 — z/zf). Certainly, when 2; 
is close to Zc the approximation of the given pump field 
fails, but the efficiency of the process should be high for 
such z. 

Fig. a shows the fields’ intensities calculated for zero 
incident first field (b) and small (but not zero) incident 
seeding second field (c); in contrast to the generation of 
two photons, here the propagation equations are not sym¬ 


metric with respect to the fields. In panel (b) one can see 
that initially the intensity of the first field grows quadrat¬ 
ically, so the solution describes difference frequency gen¬ 
eration. For the case presented in panel (c) initially the 
second field grows exponentially from the seeding. The 
following hyperbolic growth of the intensities is common 
for both cases. Again, the solutions in the transition re¬ 
gion can be less smooth depending on the initial phases 
of the fields and phases of the induced susceptibilities. 

Stimulated scattering. The schematic of the Stokes 
wave generation due to the high-order stimulated scat¬ 
tering is presented in the figure Similar to the pertur¬ 
bative case (see [1], paragraph 10.3 ), nonlinear polariza¬ 
tion can be written as: P^^(wi) = «:^^(a;i)i?i(a;i)where 
k^^(clIi) is the induced Raman susceptibility. Note that 
here we do not specify the number of laser quanta q in¬ 
volved in the process. Practically can be calculated 
exactly in the same way as it was done above, namely 
as the limit of P{uji, £i)/£i for weak £ 1 . The propa¬ 
gation equation for the Stokes wave is: +a)£’i = 

i2TrkiK^^\uJi)£i where a is linear absorption coefficient 
for frequency Wi. Its solution is £\ = £i(0) exp{g/jz —az} 
where g^, = i2'KkirS^ (loi). The main practical advan¬ 
tage of this process for effective XUV generation is the 
absence of the phase-matching problem. Our approach 
can be useful for describing XUV amplification found in 
numerical studies [II1I271. 

In conclusion, in this paper we suggest the formalism 
of the nonperturbatively induced susceptibilities which 
allows writing the propagation equations in a form sim¬ 
ilar (but not identical) to the one used in perturbative 
nonlinear optics. Under some limitations we derive the 
analytical solutions for the propagation equations, de¬ 
scribing several high-order optical properties. In particu¬ 
lar, the super-exponential (hyperbolic) growing solutions 
are found for the three-photon parametric amplification, 
never studied earlier to the best of our knowledge. The 
numerically found susceptibilities are too low to make 
these processes observable in small dilute targets like gas 
jets, but in larger targets with higher densities (high- 
pressure gas filled capillaries) such processes are feasible, 
opening a way of efficient generation of coherent XUV. 
Finally, we discuss published numerical and experimental 
results where some processes from the nonperturbative 
nonlinear optics might have already been observed. 

This study was supported by the ’’Extreme fields and 
applications” program of the Presidium of RAS. 
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Appendix I. Theory 

For completeness in this Appendix we shall start deriving equation Q from the very beginning. The microscopic 
polarization P{t) is a functional of the external field E = Eq + Ei 


P{t) = <^>[E] = $[Ao + El] 

The functional can be expanded in the Taylor series |12j (for some physical applications of such expansion, see 

[mn]): 


p{t) = p^°\t)+ + p^^\t) + ■■■ 


where 


p(o)(t) = $[F;o], 


( 5 ) 




E=Eo(t-T) 


Ei{t — r) + c.c. 




E=Eo(t-T') 


Ei(t — T)Ei{t — t') + c.c. 


here are the functional derivatives of the functional $ over the function E. 

Let us consider fields 


( 6 ) 

( 7 ) 


Eo,i = fo,i exp(—iwopt) + c.c. 


Periodicity of the field Eq allows the expansion of the functional derivatives in 
series. 


equations §, 0 in the Fourier 


5$ 

SE 


E=Eo{t-T) 


= ^G^q\^o.^G,T)eyip{-iqujii{t-T)) 


1 52 $ 
2jl^ 


= X! (fo, Wo, t') exp{-iqu;o{t - t')) 

E=Eo(t-T') q 


Substituting these expansions into equations 0,0 we have: 


p(i) (t) = J2„ £1 exp(—i(gwo + 0 Ji)t)Kq'^^'^ + £l exp(—— a;i)t)Kq + c.c. 


( 8 ) 

(9) 


( 10 ) 


P^^\t) = exp{-i[qujQ + 2 wi)t)K[+ 2 ) _|_ (£:*)2 exp(-z(ga;o - 2a;i)t)«:^ + exp(-*ga;ot)K^°’^^ + c.c. (11) 

9 


where 





dTGW(£o , Wo, r) exp(iga;o'r ± Iloit) 


( 12 ) 
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,(± 2 ) 


p+oo p+oo 

(fo, wo,a;i) = / dr dr'C^^^So , wojt') exp{iquJoT' ± iuJi{T + t')) 

Jo Jo 


(13) 


^+oo p+oo 

Kf’'^\£o,u}o,u}i)= dr dT'G^q\£o,(^o,T')exp{iqu}oT')[exp{iuJi{T-T'))+exp{-iuJi{T-T'))] (14) 
Jo Jo 

Under certain conditions (see [M. L. Pons, R. Taieb, and A. Maquet Phys. Rev. A 54, 3634 (1996)] and references 
therein) the functional 4>[Uo] in equation (H can be expanded in Fourier series according to the Floquet theorem: 


(i) = X! 4°^ (^ 0 , Wo) exp(-i9a;ot) + c.c. 


(15) 


Thus the Fourier expansion of the microscopic response P can be written as: 


P{t) = 


+ 00 +00 

^ exp{—iqujQt — imujit} + c.c. 




(16) 


q— — oo m— — oo 


where 


(fo,^l, Wo, Wl) 


+ 00 

5] «('"'^-)(£:o,wo,wi)|£-iP' 


1 = 0 , 2 . 


(17) 


(m,0) _ (m) 

Kg — Kg 

From equation ( |17| ) one can see that terms with j = 2,4,... can be understood as those describing processes where 
one, two, etc. photons from the field Ei were absorbed and the same number of photons of this field were emitted. 
Now let us consider the fields 


-Uo,i = £o,i exp(-ia;o,it + koqr + Po,i) + c.c. 


(18) 


with the real amplitudes £o,i. 

Assuming the locality of the microscopic response and substituting —wo,i^ 
we have 


—Wo,it + ko,ir + (/9o,i in equation (16) 


+ 00 +00 

Pit) = ^ ^ exp{zg(-wot + kor + (/3o) + iTO(-wit + kir + tpi)} + c.c. 

q— — oc) m— — oo 


Thus we have derived equation Q. 


(19) 


Appendix II. Some properties of the induced susceptibilities 


In this paper we consider only the centrally symmetric case (the medium is centrally symmetric and the fields are 
polarized in the same direction). 

(i) From the numerical TDSE solution we have found that while q is in the plateau region are comparable for 
different orders q. 

(ii) For wi << quiQ and in absence of resonances from equation (12) one can see that 


kW(wi)«4 


(20) 




From equation (13) for Kq 
muji << quJo 


( 2 ) 


and from similar equations for higher-order induced susceptibilities we have under 


( 21 ) 

For TO = 1, 2 this can be seen from Fig. Moreover, for higher to the difference between and Hq is visible. 
(Hi) In the following studies we shall write explicitly the frequency of the response, using the following notations: 

= quJo±uJi) = 

= qujQ ± 2uJi) = 

The permutation symmetry of the second order nonlinear susceptibilities is [T] 

^(2)*(a; = wi -I- W 2 ) = (wi = —a ;2 -I- w) = {lo 2 = uj — uji) (22) 

In the approximation of the given pump field the first equation does not have its analog for the induced suscepti¬ 
bilities. However, the second equation has the analog: 


= gwo - Wi)| = = -UJ 2 + quJo)\ (23) 

Note that the equation takes place only for the absolute values. We shall derive this equation in Appendix III. It can 
be also derived directly from the Manley-Rowe relations as it is done in [L.D. Landau, E.M. Lifshitz, Electrodynamics 
of Continuous Media (Volume 8 of A Course of Theoretical Physics), Pergamon Press, I960, §108] for y. 


Appendix III. Parametric generation of two photons. 


Let us consider the process in which q photons from the initial field are converted into one loi photon and one 
L 02 photon, see Fig. In this case q is an even number (note that for q = 2 this is a well-known process of 
four-wave mixing; however, for higher q the analogy with the four-wave mixing is hardly helpful). In the plane- 
wave approximation and slowly-varying amplitude approximation E 12 = £’ 1 ^ 2 ( 2 :) exp{i(A:i_ 2 - 2 : —Wi 2 t)}and the complex 
amplitudes of the fields are described with the following propagation equations: 


d£i,2 

dz 


x27ruj-i o T 

- -^P^-^(wi, 2 , z) exp{-z(fci_ 2 Z - Wi ,20 

A:i, 2 c 2 


where is the nonlinear part of the polarization. For to = — 1 we have from equation Q 

= K^q\wi^2 = qujQ - uj2,i) eyi^{iq{koz - ujot + (po)}E2,i 


Substituting this equation into equation (24) and assuming in the denominator fci 2 ~ ^ 1 , 2 !^ we have: 


ggl ,2 

dz 


= i2'Kki2H^q'^ (wi_2)£^2 1 exp{f Afcz -I- iqipt)} 


(24) 


(25) 


(26) 


where Afc in the considered case is Afc = qko — ki — k 2 - 


Propagation equations (26 1 are similar to those describing the parametric amplification in the usual (perturbative) 
nonlinear optics (see [T], part 9.1). Note that this similarity can also be seen from the fact that the perturbative 
nonlinear response in the presence of the pump field E 3 with the frequency 0 J 3 = wi +W 2 '- { 101 ^ 2 ) = 1 E 3 

can be converted in equation (25) substituting G'izp{iq{kQZ—ujQt+ipQ)}. However, an important difference 

from the perturbative case is that the susceptibility is complex, and its phase is not negligible. For instance, the 
phase of n^q'’ is well-understood within the recollision picture (see [ 23 ] and references therein). 

Let us introduce: 
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Equations (26) are rewritten as: 


Ai^ 2 = i27rki^2K^^\oJi,2) 


(27) 


d£i,2 

dz 


^ 1 , 2^2 1 expjzAfcz + 


Below we shall consider the case of the exact phase-matching A/c = 0. Let us introduce: 


^1,2 = Mi,2exp{zi^i^2} 


(28) 


^1,2 = au2exp{i6>i,2} 

where ui. 2 , 7 >i, 2 , ai, 2 ) 6 * 1,2 are real. Introducing 


(29) 


i)} = 6i+ qipo - ipi-(p2 


we have from equations ( 281 : 


A 0 — O 2 — 0i 


^ = aiU 2 cos(V') 

^ = 02111 cos('ip + A 0 ) 

“ 1 ^ = 0'iU2Sm(-ip) 

U 2 ^ = 0201 sin(i/! -I- A9) 

In the absence of resonances the number of quanta generated at frequencies wi and UJ 2 should be equal: 


1 d(ul) 1 djul) 
uji dz L 02 dz 


From this equation and the first pair of equations (30) we find: 


Under A6 = 0 this gives 


— cos(i/)) = — cos(i/> -I- A9) 

Wi LO2 


(31) 


Oi _ 02 
UJi UJ 2 


(32) 


Under A9j^ 0 equation (31) can be satisfied only under if) = const (for 
equations (30); for different solutions ijj can be different). Thus ^ 
equations ( 3D I we find: 


every solution of the system of linear 
- = 0 and from the second pair of 


Oi _ 02 
UJl U}2 


(33) 


This, this equation is valid irrespectively on A9. From equations (27), (29), (33) we have 
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c(l)f 




This equation presents the permutation symmetry of induced susceptibilities of the first order. Let us 
= |«: 5 ^^(a;i, 2 )|, so 01,2 = 27rA:i,2|«:g^^|- 
The solution of equation (281 is: 


^ 1 , 2 ( 2 ) = 


s^iz) =£’{^(0)exp{g=^z + i(/3f(0)} 


='^^2'^(0)exp{(g=^)*z + z:^^(0)} 


£■±(0) 


ag 2 ( 0 ) +£’i(0)exp{i(6li + q(po - j’^)} 
2 


£- 2 ±( 0 )=£-±( 0 )/a 


ip^{Q) = arg 


^ ag 2 ( 0 ) +gi (0)exp{^(6<i + gipo - 


( 0 ) =01+ qifo -ipf ( 0 ) 

where = 2Tr^/kik2\Kq^'' \ expitp^, a = yjLO\jui2-, = —A0/2, tp~ = —A8I2 + tt 
F or the phase-matched process with wi = coo/q we have 

\Re{g^)\ = 27rfco|ft:^^^|| cos(A6'/2)| 


Appendix IV. Parametric generation of three photons. 

The nonlinear polarizations for this process are: 


P^^iuJi) = = quJo -UJI- U}2) exp{iq{koZ - uj^t + lpq)}EIE2 

P^^{uj 2 ) = = gwo - 2a;i)exp{ig(fcoz - wot -f v^o)} {El)^ 


Let us denote: 

Ai = i27rfcift:^^^(a;i = gwo — wi — ^ 2 ) 


A 2 = i2’nk2K^q\uj2 = quJo — 2uJi) 
The propagation equations are written as: 


d£i 

dz 


(34) 

denote 


(35) 

(36) 

(37) 

(38) 

(39) 


Ai£i £2 expjiAfcz -I- iq+o} 


(40) 
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B£ 

= A 2 (£if exp{iAkz + iqcpoj (41) 

az 

where Ak = gkg — 2ki — k 2 - Below we shall consider the case of the exact phase-matching Ak = 0. 

Let us introduce: 


^1,2 = Mi,2exp{iv?i^2} 


^1,2 = an2exp{i6>i,2} (42) 

where ui. 2 , ^^ 1 , 2 , ai, 2 , ^ 1,2 are real. Introducing 


-ip = 01 + q(po - 2(pi - ip2 


A9 = 02 — 01 


we have from equations (401, (41): 


^ = aiUiU 2 cos('0) 

^ = 02uf cosipj + A0) 

= aiUiU 2 sin('0) 
U2% =a2ulsm{-iP + A0) 


(43) 


In the absence of resonances the number of quanta generated at the frequency wi should be twice that of quanta 
at the frequency UJ 2 '- 


1 d{ul) 2 d{u2) 
uji dz L 02 dz 


From this equation and the first pair of equations (43) we find: 


— cos{ip) = 2— cos('0 -|- A0) 

Wi UJ2 


Case 1: AO ^ 0. Equation (44) can be satisfied only under ip = const: 


Thus 


tan(^/>) 


cos(A0) 


2 \K^q^ (W 2 )| 


1 

sin(A6*) 


(44) 


(45) 


dip _ 

dz dz dz 


and from this equation and the second pair of equations (43) we find: 


(46) 


ui{z)/u 2 {z) = const 

The solution of the first pair of equations (43) satisfying is: 


(47) 
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Ul,2iz) = Ui^2{0)z - r 

J. Z j 2 . 


(48) 


where 


1 


-Zc = 


U 2 ( 0 )ai cos(V') 


(49) 


and 


Ui(z) / Oi cos('0) 
U 2 (z) y 02 COs(V’ + A0) 


Taking into account equation (44), this equation can be rewritten as 


uijz) 

U2iz) 


2iL0i 

UJ 2 


(50) 


Thus only if the initial conditions satisfy this equation, the solution satisfies equation (47). Finally, from equa¬ 
tions (46) , (50), and the second pair of equations (43) we find: 


cos(A0) = 


K^q \oji)\ 2|s:^^ -'(a;2)| 

(^i)l 


( 2 )/ 


,( 2 ) 


( 0 ^ 2)1 


.( 2 ) 


(51) 


Thus we find the connection between phases and absolute values of the second-order induced susceptibilities. 
Case 2: Ad = 0. Equation (44) is satisfied for every ip] it is satisfied under: 


0-1 _ 2 


Wi 


0 J 2 


(52) 


So 


|«f(a;i)|=2|«f(a;2)| 


(53) 


This equation presents the permutation symmetry of the induced susceptibilities of the second order for Ad = 0. 


Note that equation (51) turns to (53) when Ad = 0. 
Let us denote: 


£ 1 , 2 ( 2 ) = £' 1 , 2 ( 2 ) exp{z(/3i, 2 ( 0 )} 


so that £' 1 . 2 ( 0 ) are real. Then equations 


(40), (41) can be written as: 


d£'i 

dz 


ai£'i £'2 exp{ii/>(0)} 


(54) 


= 02 (£'i)^exp{z'(/)(0)} 


(55) 


where ^^(O) = di -|- qipo - 2 vji(0 ) - (/J 2 ( 0 ) 
(/3i(0) and <^ 2 ( 0 ) are chosen so that 'ip(O) 


(note that di = 62 because we consider the case Ad 
= 0, the solution of equations (54), (55) remains real. 


= 0). If initial phases 
Substituting 


£' 1 ( 2 ) 


^^(£' 2 ( 2 ) 2 -£ 2 ) 
V W2 


(56) 
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where £ = 


equation (52)): 


^£' 2 {^Y — ^£'i(0)^, system of equations (54), (55) is presented as a single equation (note also 


d £'2 

dz 


ai{£'l-£^) 


Such equation is known as the Riccati equation. Its solution is: 


(57) 


, ~ fiexp{-gz} + expjgz} 

^^exp{-gz} - exp{gz} 


(58) 


where g = ai£, 


£'2(0)+g 
£' 2 ( 0 )-£" 


Thus equations (56) and (58) 
certain propagation distance: 


give the solution. 


The typical feature of this solution is a hyperbolic growth after a 


£'n2(^) oc 


1 

1 - z/Zc 


The solution goes to inhnity at 


_ ln(M) 

Z (2 - _ 

2 g 


(59) 


Thus we can see that in both cases A9 7 ^ 0 and A6 = 0 there is a hyperbolic growth of the solution after a certain 
propagation distance. 













